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ABSTRACT 

The growth of a turbulent spot in a laminar 
boundary layer, aa the spot evolves from a localized 
disturbance in the layer, is simulated numerically 
using a three-dimensional vortex fllamenc description 
of the vortlcitv field. The filaments are marked with 
a sequence of node points which are tracked in a 
Lagranglan reference frame. Velocity computation is 
done by Blot-Savart integration. Although some dis- 
crepancies with experiment appear to exist in the near 
wall region, the gross properties of the spot, includ- 
ing the velocities of the leading and trailing edges 
and the velocity perturbations away from the wall, are 
in good agreement with experiment. 

INTRODUCTION 

Studies of the turbulent spot in a laminar bound- 
ary layer date back to the theoretical and experimen- 
tal work of Emmons (1) and Mltchner (2), who were con- 
cerned with the role of spots in the transition pro- 
cess. Such spots are now recognized as interesting 
flow phenomena in their own right and, in addition, 
spots or spotlike objects have been suggested as the 
coherent flow structure responsible for much of the 
turbulent transport in a fully developed turbulent 
boundary layer (3,4). Recent experiments by U-^nanski, 
et al. (5) and Cantwell, et al. (6) confirmed ear'.ls. 
observations of Schubauer and Klebanoff (7) and pro- 
vided much new information on the structure of turbu- 
lent spots. In particular, entrainment rates in the 
plane of symmetry and the mean flow field within the 
spot were obtained by ensemble averaging over a large 
number of spots. In another recent experiment. Caster 
and Grant (8) studied the evolution of a wave packet 
in a laminar boundary layer produced by a localized 
perturbation at the wall boundary. The amplitude of 
the perturbation was very small so that the linear 
development of the resultant disturbance could be 
observed for some time after its initiation. Indeed, 
model calculations b> Caster (9) based on linear theory 
agree well with the e.-perlmental measurements of Caster 
and Crant (8) during toe early stages, but some quali- 
tative differences were noted in the later stages of 
growth, presumably due to nonlinear mechanisms. 

Despite these significant advances in our under- 
standing of turbulent spots, a number of questions 
remain, such as the interaction of two or more spots, 
the formation of sublayer streaks, the three-dimensional 
vortlclty field within the spot, the nonlinear mechan- 
lce of growth, and che birth of new spots. 

In this paper, the growth of a turbulent spot in 
a laminar boundary layer, as the spot evolves from a 
localized disturbance in the layer, is simulated 


numerically using a vortex filament description of the 
vortlclty field. The attempt is to provide additional 
insight into the structure and growth of a turbulent 
spot. In many respects, the method employed (10,11) is 
an extension to three dimensions of the two-dimensional 
vortex method (12,13) and is done in the spirit of the 
large-eddy simulation technique (14) in which small- 
scale turbulence is not simulated directly but la 
modeled. In an earlier numerical study Hama (15), 
motivated by his experiisental obeervatii s of boundary- 
layer transition, investigated the motion of a single 
vortex loop. 

Previous three-dimensional numerical simulations 
of wall-bonded shear flows have Included channel flow 
simulations by Deardorff (16) and Schumann (17) using 
an artificial law-of-tne-wall boundary condition and 
the recent Improved simulation by Mo in et al. (18) 
using the no-slip condition. In addition, transition 
phenomena have been simulated numerically by Orszag 
(19) and by Wray et al. 120). 

NUMERICAL METHOD 

For numerical purposes, a turbulent spot is viewed 
as a localized region of boundary-layer vortlclty that 
has undergone three-dimensional displacements, strains, 
and rotations from its original monodlrectlonal laminar 
state (see Fig. 1). Each tube of vortlclty in the 
boundary layer is decomposed into its straight, unper- 
turbed contribution and into a loop representing its 
contribution to the perturbed vortlclty field, as shown 
in Fig. 2. The velocity field is decomposed, accord- 
ingly. as 

u(x) - U(y)e x + u* (x) 

where U(y) is the laminar profile and u' is given by 
a sum of Blot-Savart integrations over all the filament 
curves, C^, plus their images. The image contributions 
ensure tangency of the flow at the wall. In this ini- 
tial study we ignore the generation of vortlclty at the 
wall boundary due to the disturbance, i.e., we use the 
lnvlscld boundary condition at the wall. Blot-Savart 
integration yields an exact representation of the bound 
ary conditions at infinity in all directions. 

Each filament is represented by a space curve 
X|(t,t), where C la a parameter along the curve; by a 
circulation T i ; and by an effective core radius cj, 
which parameterizes the assumed Gaussian vortlclty dis- 
tribution within the filament. Thue the perturbed vor- 
tlclty field is given by 

w(x,t) • Tj f exp[-|x-x 1 | 2 /o 1 J ]/(/so 1 ) 1 -^- dt 


*Presented at the Second Symposium on Turbulent Shear Flows, July 2-4, 
1979, Imperial College, London, England. 
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A filament or curve represents a collection or vortex 
line* to that the filament velocity must be computed 
as an appropriate average over the filament core. We 
use the following equation of notion for the space 
curvet: 



♦ # X U (y ^ ) ♦ Image contributions 

The equation is seen to be the Blot-Savart induction 
law for an Infinitesimal filament but corrected In the 
near field for finite core size. Using a • 0.2065 
we recover the correct speed of a ring vortex in the 
limit o i << ring radius. 

For numerical purposes, the fi lamer ca are marked 
with a sequence of node points which are tracked in a 
Lagrangian reference frame. The space curve x^(C,t) 
la assumed piecewise linear, passing through each node 
point. The spacing between node points Is O(o^) to 
minimize discretization error. For example, satisfy- 
ing this constraint yields a computed speed of a ring 
vortex that is essentially independent of the number 
of node points. 

Discretization in the directions normal to the 
vorticity vector oust also be treated carefully. If 
the nominal spacing between vortices is so large that 
adjacent cores do not overlap, the smooth vorticity 
field that exists in the undisturbed regions near the 
spot Is poorly approximated. Given a number of hori- 
zontal layers of vortex filaments, equally spaced in 
the vertical direction, an optimum core size as well 
as s set of optimum circulations for each horizontal 
layer is determined that minimizes the mean square 
error in the desired vorticity distribution. The 
result of this procedure is displayed in Fig. 3(a) for 
the esse of eight layers, whereas the effect of choos- 
ing a smaller value of core size is shown in Fig. 3(b). 

For many incompressible flows in which the vor- 
ticity Is concentrated in small regions, vortex methods 
have the conceptual advantage that coraput at tonal points 
need only appear in the vort icit v-containing rev ions. 

In thia problem only the region of perturbed vortlnty 
requires points. Thus, in this simulation, points are 
added automatically In the atreamvlse and. In some 
cases, the spanwlse direction, as the spot grows. More 
specifically, the maximum perturbation of x* In the 
farthest upstreau and farthest downstream filaments in 
each layer is monitored. As soon as that perturbation 
reaches a critical value (usually 0.004 of the dis- 
placement thickness, 6*) a new, adjacent filament is 
added. This feature, which helps minimize the number 
of computational points, Is particularly Important 
because O(N^) operations per time step are required 
to confute the velocity at each of N vortex node 
points. For reasonable total run times (1-2 epu hr), 
the maximum number of points representing one-half the 
spot la <2000 on the CDC 7600 and <8000 on the 
IUUc IV. 

The range of computed scales of motion is quite 
large; e.g., in the streamvlse direction the ratio of 
the largest to the smallest computed scale is ^125 to 
1. Nevertheless, this range is not large enough to 
span those scales found experimentally (see, e.g., 
aluminum particle photographs in Ref. 6). Thus, the 
effect of the smaller or subgrid scales on the com- 
puted scales must be modeled. The usual procedure 
used in finite-difference simulations is to use a sub- 
grid viscosity dependent on the local strain rate (14). 


In this study we simply maintain a constant filament 
core size. Indeed, this procedure is basically equiva- 
lent to using a subgrid viscosity dependent on the 
local rate of vortex stretching. 

RESULTS 

Figure 4 displays the computed grid of vortex 
filaments as the spot develops. All the results given 
in this paper were obtained on the Illiac IV, although 
development of the algorithm was done on a CDC 7600. 

As shown at t • 0, the initial perturbation imposed on 
the boundary layer consists of a localized deformation 
of the vortex lines in the streamwlse direction with a 
maximum displacement equal to 0.3. The net effect of 
this dtformarlon is equivalent to adding a ringlike 
blob o» vorticity to the laminar boundary layer, 
or lent et so that the blob la initially moving upward 
away frot the plate. The maximum streamwlse velocity 
perturbation in this initial disturbance is -0.07 U M . 
Velocities and lengths are normalized so that U. - 1 
and 3U/3y| wa n ■ 1. The unperturbed vorticity profile 
is that shown in Fig. 3(a) and the corresponding 6* 
is 1.21. Eight levels of filaments were used in the 
v-direction normal to the piste and were spaced 0.2 
apart. The same spacing was used in the streamwlse (x) 
direction due to the isotropy of the velocity field 
induced by each filament. The Initial syaclng of points 
along each filament was 0.67. As the simulation pro- 
ceeded, however, points were automatically redistrib- 
uted to provide a higher density of points in regions 
of high curvature. The lifting of vortex lines is evi- 
dent in the side view of Fig. 4. The stretching of 
vortex lines and the development of streamwlse vortic- 
ity is quite evident in the plan view. In the simula- 
tion shown in Fig. 4, 2184 points in 104 filaments ware 
used at t • 0. At the end of the run of 71 steps 
(At - 0.3, first-order Euler explicit time differenc- 
ing), there were 12,222 points in 582 filaments repre- 
senting the spot. The spot was assumed to be symmetric 
with respect to the z • 0 plane. Therefore, the 
number of computational points requiring nodal veloci- 
ties was only approximately half the total number 
representing the disturbance. 

In Fig. 5, a time sequence of contour plots Is 
shown for the perturbed streamwlse velocity in the 
plane of symmetry. The plots are dominated by a growing 
patch of streamwlse velocity deficit with a peak of 
-0.37 U,. at the final time, t • 21. Small regions of 
positive perturbation are evident near the wall in the 
later »ges. Corresponding plots of vertical velocity 
and perturbed spanwlse vorticity are shown in Figa. 6 
and 7, respectively. As shown, the vertical velocltlea 
are predominately positive. The vorticity plot shows 
an elongated region of negative vorticity riding on a 
similar region of positive vorticity. 

To compare with experiment, profiles of the spot 
boundary were determined from plots similar to those in 
Fig. 5 according to the criterion of Ref. 5: 2% devia- 

tion of the streamwlse velocity from its free-stream 
value. Shown in Fig. 8 are trajectories of the leading 
and trailing edges measured In the plane of symmetry et 
y • 0.6, 0,9, and 1.2. (These y values have undis- 
turbed laminar velocities of 0.55 U„, 0.75 U w , and 
0.89 !’„, respectively.) The trajectories for y - 1.2 
yield leading- and trail lng-edge velocities of 0.88 U. 
and 0.48 U -f respectively. These values compare well 
with those of Wygnanskl et al. (5), who found 
^LE * 0.89 U„ and Uje • 0.50 U*, and with those of 
Schubauer and Klebanoff (7), who reported Ule» 0.88 U. 
and Ute " 0.50 U„. Tral 1 ing-edge velocities, however, 
show a significant variation with y. The plot also 
reveals the virtual origin (x 0 , t 0 ) of the epot. 
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Velocity histories at fixed point* arc known In 
Fig*. 9 ana 10. The u-veloclty record* aw-.y from the 
wall show the cherecterl*t lc dip end recovery In agree- 
ment with the measurement* of Kefs. 5 and t> ’used on 
ensemble mean* over large number* of spot). The simu- 
lated u-velocltie* In the plane of symmetry near the 
wall show no Indication of the hlgh-.,peei rise which 1* 
observed experimentally. The normal velocity record 
shown In Fig. 10 has an Initial dip followed by a rise 
to positive values in agreement with the ensemble- 
averaged tesults of Ref. 5. 

Streamwlse vciocltv pert ut oat ions In three hori- 
zontal plane* at t • 18 art shown In Fig. 11. 

Although the velocity histories taken In the plane of 
Nyasoet ry (Fig. 9) showed no evidence of htgn-spced 
fluid near the wall, the contour plot for y • 0.3 
■how* extensive region* of positive perturbation away 
frosi the plane of symmetry. 

CONCLUSIONS 

The numerical simulation of a spotlike disturbance 
In a laminar boundary layer using vortex filaments has 
provided Insight Into the development of the three- 
dimensional vortlclty and velocity fields within the 
apot. Cross properties of the spot away from the wall. 
Including the velocities of the leading and trailing 
edges and the velocity perturbations, are in good 
agreeswnt with experiment. 

Discrepancies with experiment appear to exist In 
the near wall region, especially at the trailing edge. 
This suggests that future simulations should include 
the efface of vortlclty production at the wall and an 
Improved representation of viscous dissipation. 

Another source of disagreement may be the relatively 
shore history of the computed spot. Experimental 
results (5-7) are based on spots that have traveled 
several hundred 4* whereas the computed spot has 
traveled about 20 6*. Another possibility Is that the 
experimental data were obtained from ensemble averaging 
over a large number of spots, a procedure that possibly 
filters out Important coherent nx>t Ions on the scale of 
several 4* or more In addition to unwanted "noise." 
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Figure 1.- Spot geometry. 
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Figure 2.- Decomposition of vorticity field. 
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(a) Optimum vortex representation. 



y 


(b) Suboptimum choice of core size. 
Figure 3.- Vorticity distributions. 
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Figure 6.- v' contour plot; same contour levels as figure 5; t 


contour plot; contours 0.1, 0.3, 0.5 
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Figure 8.- Trajectories of leading and trailing edges 
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Figure 9.- u vs time, x = 8.0, z = 0.0. 


9 





Figure 11.- u' contour plots in horizontal planes; 
contour levels ±2%, ±6%, ±10%, . . . 
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